A new formula for the saddle-to-scission time

We propose a new formula for the saddle-to-scission time that is more general that the one based on Kramers' approach. Its validity and applicability is then studied in detail. Such a formula is useful for the evaluation of the fission time of very heavy nuclei.     

Assessment of some iterative methods for non-symmetric linear systems arising in computational fluid dynamics

Various tests have been carried out in order to compare the performances of several methods used to solve the non-symmetric linear systems of equations arising from implicit discretizations of CFD problems, namely the scalar advection-diffusion equation and the compressible Euler equations. The iterative schemes under consideration belong to three families of algorithms: relaxation (Jacobi and Gauss-Seidel), gradient and Newton methods. Two gradient methods have been selected: a Krylov subspace iteration method (GMRES) and a non-symmetric extension of the conjugate gradient method (CGS). Finally, a quasi-Newton method has also been considered (Broyden). The aim of this paper is to provide indications of which appears to be the most adequate method according to the particular circumstances as well as to discuss the implementation aspects of each scheme.     

Incremental unknowns preconditioning for solving the Helmholtz equation

An efficient preconditioner is developed for solving the Helmholtz problem in both high and low frequency (wavenumber) regimes. The preconditioner is based on hierarchical unknowns on nested grids, known as incremental unknowns (IU). The motivation for the IU preconditioner is provided by an eigenvalue analysis of a simplified Helmholtz problem. The performance of our preconditioner is tested on the iterative solution of two‐dimensional electromagnetic scattering problems. When compared with other well‐known methods, our technique is shown to often provide a better numerical efficacy and, most importantly, to be more robust. Moreover, for the best performance, the number of IU levels used in the preconditioner should be designed for the coarsest grid to have roughly two points per linear wavelength. This result is consistent with the conventional sampling criteria for wave phenomena in contrast with existing IU applications for solving the Laplace/Poisson problem, where the coarsest grid comprises just one interior point. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Solving large-scale semidefinite programs in parallel

We describe an approach to the parallel and distributed solution of large-scale, block structured semidefinite programs using the spectral bundle method. Various elements of this approach (such as data distribution, an implicitly restarted Lanczos method tailored to handle block diagonal structure, a mixed polyhedral-semidefinite subdifferential model, and other aspects related to parallelism) are combined in an implementation called LAMBDA, which delivers faster solution times than previously possible, and acceptable parallel scalability on sufficiently large problems. This work was supported in part by NSF grants DMS-0215373 and DMS-0238008.     

Identification of a dissipation coefficient by a variational method

A generalized inverse problem for the identification of the absorption coefficient for a hyperbolic system is considered. The well-posedness of the problem is examined. It is proved that the regular part of the solution is an L ₂ function, which reduces the inverse problem to minimizing the error functional. The gradient of the functional is determined in explicit form from the adjoint problem, and approximate formulas for its calculation are derived. A regularization algorithm for the solution of the inverse problem is considered. Numerical results obtained for various excitation sources are displayed.     

STABLE SOLUTION OF TIME DOMAIN INTEGRAL EQUATION METHODS USING QUADRATIC B-SPLINE TEMPORAL BASIS FUNCTIONS

This paper is concerned with stable solutions of time domain integral equation （TDIE） methods for transient scattering problems with 3D conducting objects. We use the quadratic B-spline function as temporal basis functions, which permits both the induced currents and induced charges to be properly approximated in terms of completeness. Because the B-spline function has the least support width among all polynomial basis functions of the same order, the resulting system matrices seem to be the sparsest. The TDIE formula-tions using induced electric polarizations as unknown function are adopted and justified. Numerical results demonstrate that the proposed approach is accurate and efficient, and no late-time instability is observed.     

一种新的高温超导磁体电磁结构优化设计方法

本文阐述了用于高温超导磁体电磁结构优化设计的一种新方法,与一般优化方法不同之处在于其不用定义磁体的初始结构形状,采用该方法可以得到满足规格和性能要求的较佳磁体结构.通过此方法优化设计了一个中心磁场为3T的Bi系高温超导磁体,其设计结果与满足同样要求且经过优化的单螺管高温超导磁体的设计结果作了对比,由结果表明该方法具有很好的优化效果.     

On the convergence rate and optimization of a numerical method with splitting of boundary conditions for the stokes system in a spherical layer in the axisymmetric case: Modification for thick layers

The convergence rate of a fast-converging second-order accurate iterative method with splitting of boundary conditions constructed by the authors for solving an axisymmetric Dirichlet boundary value problem for the Stokes system in a spherical gap is studied numerically. For R/r exceeding about 30, where r and R are the radii of the inner and outer boundary spheres, it is established that the convergence rate of the method is lower (and considerably lower for large R/r) than the convergence rate of its differential version. For this reason, a really simpler, more slowly converging modification of the original method is constructed on the differential level and a finite-element implementation of this modification is built. Numerical experiments have revealed that this modification has the same convergence rate as its differential counterpart for R/r of up to 5 × 10³. When the multigrid method is used to solve the split and auxiliary boundary value problems arising at iterations, the modification is more efficient than the original method starting from R/r ～ 30 and is considerably more efficient for large values of R/r. It is also established that the convergence rates of both methods depend little on the stretching coefficient η of circularly rectangular mesh cells in a range of η that is well sufficient for effective use of the multigrid method for arbitrary values of R/r smaller than ～ 5 × 10³.     

General Aspects and Optimization of Enantioselective Biocatalysis in Organic Solvents: The Use of Lipases [New Synthetic Methods (76)]

Enantioselective biocatalysis in nonaqueous media is becoming increasingly important in preparative synthetic chemistry. This article discusses (1) the general catalytic properties of enzymes in nonaqueous environments, (2) the basic principles that govern lipase-catalyzed enantioselective esterification and transesterification reactions in organic media for the preparation of optically active acids and alcohols, (3) the determination of kinetic and thermodynamic parameters, and (4) the quantitative analysis of published data.     

Homogeneous immunoassays

Applied Biochemistry and Biotechnology -